Strongly interacting photons in hollow-core waveguides 



O 
(N 

o 

Q 



Ephraim Shahmoon,^ Gershon Kurizki/ Michael Fleischhauer,^ and David Petrosyan^''^ 

'Department of Chemical Physics, Weizmann Institute of Science, Rehovot 76100, Israel 
^Department of Physics and Research Center OPTIMAS, 
Techmsche Universitdt Kaiserslautern, D-67663 Kaiserslautern, Germany 
'^Institute of Electronic Structure & Laser, FORTH, GR-71110 Heraklion, Crete, Greece 

(Dated: December 17, 2010) 

Hollow-core photonic-crystal waveguides filled with cold atoms can support giant optical non- 
linearities through nondispersive propagation of light tightly confined in the transverse direction. 
Here we explore electromagnetically induced transparency is such structures, considering a pair of 
counter-propagating weak quantum fields in the medium of coherently driven atoms in the ladder 
configuration. Strong dipole-dipole interactions between optically excited, polarized Rydberg states 
of the atoms translate into a large dispersive interaction between the two fields. This can be used to 
attain a spatially-homogeneous conditional phase shift of n for two single-photon pulses, realizing a 
deterministic photonic phase gate, or to implement a quantum nondemolition measurement of the 
photon number in the signal pulse by a coherent probe, thereby achieving a heralded source of single 
or few photon pulses. 

PACS numbers: 42.50.Gy, 42.65.-k, 03.67.Lx, 32.80.Ee, 
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Photons are ideal carriers of information in terms of 
transfer rates and distances. Yet, scalable and efficient 
quantum information processing [l| with photons would 
require implementing deterministic quantum logic be- 
tween single-photon qubits Q , which is hindered by the 
weakness of optical nonlinearities in conventional media. 
Highly enhanced nonlinear interactions in atomic vapors 
Q in the regime of electromagnetically induced trans- 
parency (EIT) [3-0] have emerged as a promising route 
to circumvent these difficulties and to achieve large condi- 
tional phase shifts <f> for pairs of slowly propagating pho- 
tons. Attaining the phase shift of = tt would amount to 
realizing the universal CPHASE gate for photonic qubits 

i- 

Among the many relevant proposals 0-0 , one of the 
most promising schemes is based on employing EIT in 
a ladder configuration [ll'], wherein the photon-photon 
interaction is mediated by strong dipole-dipole interac- 
tions (DDIs) between optically excited Rydberg states 
of the atoms [3, [l^. An important advantage of this 
scheme is that the long-range nature of the DDI relaxes 
the need for tight focusing of the quantum fields to the 
atomic absorption cross-section <j ~ , which is close to 
the diffraction limit. 

In Ref. [ll| we have presented an effective one- 
dimensional (ID) treatment of the dynamics of 
two slowly counter-propagating, weakly-focused single- 
photon pulses. We have done so by considering the elec- 
tric fields only on the propagation axis, and have shown 
that, for a pair of photons passing through each other, 
the accumulated conditional phase shift (f) can be both 
large and uniform in the longitudinal direction. In free 
space, however, the ID treatment of interacting quantum 
fields is incomplete as it does not capture the diffraction 
effects and the fact that, in the transverse direction, the 
resulting phase shift is inhomogeneous due to the rela- 
tive coordinate dependence of the DDI potential [1^1 ■ To 



remedy these problems and achieve non-diffracting, uni- 
form transverse phase-fronts, here we propose to impose 
onto the quantum fields only a single transverse mode by 
confining them into a hollow-core photonic-crystal waveg- 
uide [iTlflsj filled with an ensemble of cold alkali atoms 
[l^ . In what follows, we present a rigorous derivation 
of ID propagation equations for two interacting quan- 
tum fields. We extend our earlier scheme by considering 
the atomic level configuration involving different Ryd- 
berg states. We then discuss the conditional phase shift 
for two single-photon pulses. Furthermore, we analyze a 
quantum nondemolition measurement (QND) of the pho- 
ton number in the signal pulse inducing a phase shift of 
the coherent probe pulse. This can serve as a heralded 
source of single or few photon pulses. 

We begin by assuming that the transverse intensity 
profile of the counter-propagating fields Ei and E2 in 
the cylindrically symmetric waveguide is described by 
a Gaussian e~^^/^s of width Wf, where r±_ = \r±\ is 
the distance from the propagation z axis. The corre- 
sponding electric field can then be expressed as Ei{r) = 
eie-''±/^'"}£i{z) {I = 1,2), where ei = y^huji/2eoV is 
the field per photon of frequency uji within the quanti- 
zation volume V = nw'jL, with L the waveguide length, 

while £i{z) = J2 
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is the traveling-wave field op- 
erator, given by a superposition of bosonic operators af 
for the longitudinal field modes k, yielding the commu- 
tation relations [£i{z), £j,{z')] = Ldu'5{z — z'). An en- 
semble of N cold atoms is trapped in the hollow core of 
the waveguide tl£|; the corresponding atomic density is 



then p{y) = (ttw^) 



LA 



''{N/L), where Wa{^ Wf) is 



the width of the transverse Gaussian distribution. The 
level configuration of the atoms, all of which are initially 
prepared in the ground state \g), is schematically shown 
in Fig. [ija). The quantum fields i?i.2 resonantly inter- 
act with the atoms on the transitions |g) — > |ei^2), re- 
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FIG. 1: (a) Level scheme of atoms resonantly interacting with 
quantum fields i?i,2 and classical driving fields r2i,2 on the cor- 
responding transitions. Vdd denotes the DDI between atoms 
in Rydberg states \d). (b) The quantum fields transversely 
confined in a hollow-core waveguide of length L filled with the 
atoms, counterpropagate as dark-state polaritons ^1,2 having 
slow group velocities ui,2 and interacting via long-range po- 
tential A12 (21—2:2) mediated by Vdd. (c) The potential Ajj' (C) 
of Eq. ((51), as a function of dimensionless distance C,, in units 
ol2Cu'IK^wf Hz. 



spectively. The intermediate states |ei_2) are resonantly 
coupled by two strong (classical) driving fields with Rabi 
frequencies ili_2 to the Rydberg states |(ii,2)- In a static 
electric field E'st^z, these Rydberg states possess perma- 
nent dipole moments p — Inqeage^, where n and q are 
the (effective) principal and parabolic quantum numbers, 
e is the electron charge, and oq is the Bohr radius (2^. A 
pair of atoms at positions r and r' excited to states \di) 
and \dii) interact with each other via a DDI potential 
Vdd resulting in an energy shift 
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where ■& is the angle between vectors and r' — r, and 
Civ = pd, prf,, /(47reo) is proportional to the product of 
atomic dipole moments pd, = p jd;). We assume that 
state mixing within the same n manifold is suppressed by 
a proper choice of parabolic q and magnetic m quantum 
numbers [20j . 

We use collective atomic transition operators ct^^ (r) = 
l/-^rX]j=i averaged over the volume element 

AV containing Nr = p{r) AV » 1 atoms around posi- 
tion r. In the frame rotating with the frequencies of the 
optical fields, the interaction Hamiltonian H = V"af + Vdd 
contains the atom-field and DDI terms 

T4f = -hfd\p{r)Y, [ffie-'--/2-?£K2)'^e,<,(r) 
1=1,2 

-Kf7/CTd,e,(r)] +H. c, (2a) 

Fdd = njd\p{Y) jd\'p{v') 

E '^d,d,(r)An'(r-r')<Td,,d,,(r'), (2b) 

i,i' = l,2 

where gi = {pgeJTi)ei is the corresponding atom-field 



coupling constant, with p^g, being the dipole matrix el- 
ement on the transition |g) — >■ |e;). 

Using Hamiltonian H, we derive the Heisenberg- 
Langevin equations for the atomic operators (Tgej(r), 
dgdi (r) and the propagation equations for the slowly- 
varying quantum fields £i{z). Solving for the atomic op- 
erators perturbatively in the small parameters giSi/fli 
and in the adiabatic approximation 0, [^, [ll| , and after 
substituting into the equations for the fields, we obtain 
the following propagation equations for the dark-state 
polaritons = c/vi £1 t^]. 



-isin^eiSi{z,t)-^i{z,t), (3) 



the sign or "— " corresponding to / = 1 or 2, respec- 
tively, vi = c cos^ 01 is the group velocity of the corre- 
sponding field in the EIT medium, and the mixing angles 
9i are defined through tan^ 6*; = (gf A^/|r2;p)(w/w^)2^ 
with ui = WaWf{w'^ + w'^j)^^^'^ . Operators Si{z,t) are 
responsible for the self- and cross-phase modulation be- 
tween the fields, 

Si{z,t) = y [ dz'[Aa{z- z')smHiii{z',t) 
^ Jo 

+Au'{z-z')sm^eidi,{z',t)], (4) 

where 2/ = '^j'^i = {c/vi)£j£i are the polariton intensity 
(excitation number) operators in the EIT medium, which 
correspond to the photon number operators outside the 
medium (w; = c) [4], while the effective one-dimensional 
DDI potentials Aw (z — z') result from Aw (r — r') upon 
double integration over the transverse coordinates, 



Aw{z-z') 
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(ttui^)^ 



d'rxyd'r^e-('-i+'-")/'^ A,,,(r - r') 
[2|C|-y^(l + 2CV'erfc(|C|)],(5) 



C= {z-z')/V2w. 

As seen in Fig. [TJc), Aiii{Q is sharply peaked around 
C = with the range (FWHM) of - 0.65. 

It follows from Eq. ([3|) that the intensity operators Ii 
are constants of motion: Ii{z, t) — Ii{z^vit, 0), the upper 
(lower) sign corresponding to I = I {I = 2). The solution 
for the field operators then reads 



i'i{z,t) = exp 



-i sin^ 01 / dt'Si{z^vi{t-t'), t') 
Jo 



i'i{zTvit,0). 
(6) 



The validity of this dissipation-free solution hinges on the 
following assumptions: (i) The duration T; of each pulse 
exceeds the inverse of the corresponding EIT bandwidth 
Suji = \ / {'jgei\/ KiL), where 7gei is the transversal re- 
laxation rate and m Qp is the resonant absorption co- 
efficient, with Q = 3A^/(27r) the absorption cross section 
on the transition |.g) — >■ |e;) and p ~ iV/[7r(w^H-?i;^)L] the 
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X {n2\ exp 



dt' I dz'Ai2{z- z' ~v{t~-t')) 
n2). (10) 



"'0 

xi2{z' + vt',0) 



effective atomic density. With vi = 2|rj;p/(K/7geJ, this the expectation value of the probe field we have 
yields the condition {kiL)~^/^ ^ Tivi/L < 1 which re- 
quires a medium with large optical depth niL ^ 1 [1, [^. 
(ii) The DDI induced frequency shifts lie within the EIT 
bandwidths, sin^ 6'z(S',(z)) < <5wz, Vz S [0,L]. (iii) The 
propagation/interaction time of each pulse tout — L/vi is 
limited by the relaxation rate 7^^; of the iTgd; coherence 

via toutjgd, < 1- 

In what follows, we employ Eq. ^ to demonstrate the 
quantum phase gate between two single-photon pulses 
-£'1,2, and to realize a quantum nondemolition measure- 
ment of photon number in the signal pulse E2 by a co- 
herent probe pulse Ei . For simplicity of notation, we set 
9^,2 = 9, i.e.,gfN/\ni\^^giN/\n2\'. 

We are concerned with the evolution of input state 
l^in) = |li) II2) composed of two single-excitation 
wavepackets |1;) = [-^ /dz//(z)^|(z)] |0) whose spa- 
tial envelopes inside the medium fi{z) — (0| 5';(z,0) \ 



As before, we assume that at t = the probe and signal 
pulses are centered, respectively, at z = and z = L. 
The output probe field at tout — L/v and z — L is then 



-ai(0) 



X (7121 exp 



sin^e 



L/v pL 

dt' dz'Ai2{z' -vt') 



are normalized as Jdz\fi{z)\'^ — 1. With the opera- 
tor solution (O, for the (equal-time) correlation ampli- 
tude or the "two-photon wavefunction" Fi2{zi, Z2,t) = 
(0| *i(zi,i)*2(^2,i) |$in) we obtain 

-£"12(2:1,^:2,^) = fl{zi -Vt)f2{z2 + Vt) exp[i(j)i2{zi, Z2,t)],{7) 

(t)i2{zi,Z2,t) = -sm'^9 [ dt'Ai2{zi- Z2-2v{t-t')). (8) 



^0 

xi2{z' + vt',0) 



\n2). (11) 



Hence, the two polaritons counterpropagate in a shape- 
preserving manner with group velocities ±v. Since 
Xi'^i \ = 0, the self-interaction within each pulse is 
absent, while the cross-interaction between the pulses re- 
sults in the phase-shift ([8]) . Assume that at t = the first 
pulse is centered at zi = and the second pulse at Z2 = 
L, while after the interaction, tout = L/v, the coordinates 
of the two pulses are zi = L and Z2 = 0, respectively. 
The accumulated phase-shift is then (/)i2(i, 0, i/u) = 

— siii^9/v Jf^dz' A12 {2 z' — L). To evaluate the integral, 
we replace the variable (2z' — L)l\/2w — > C' and extend 
the integration limits to Lji^splvS) — 00, obtaining 



:^12 



C12 sin^ 9 



(9) 



which is spatially uniform and the state of the system 
at tout is |<&out) = e"^" \^in)- Since for input states 
1^2) (™, n = 0,1) there is no phase shift when m + 
n < 2, the conditional two-photon phase shift (f>i2 = tt is 
equivalent to the cphase gate |<I>out) = (-1)™" \mi) \n2) 



equ] 
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We next consider the probe pulse in a multimode co- 
herent state |ai) = Hk \cti), which is an eigenstate of the 
field operator '^i{z) with eigenvalue ai{z) = a^e'''^: 
\E'i(z) |q;i) = ai{z) \ai). The signal pulse can be in any 
superposition or mixture of the n-photon number states 
|n2) = J dz f2{z)^l{z)y^ |0). Given an input state 

l^in) = |n2), and neglecting the self-interaction, for 



Recall that the DDI potential Aw (z) is sharply peaked 
around z = with the range <5z < w <C i [Fig. [U^c)], 
while J^^dzAiii{z) = —2Cw/{1iw'^). On the other hand, 

in the EIT medium, (n/|I;(z) \ni) = rt|/;(z)p are smooth 
pulses of length Tiv ^ L. To evaluate the integral in the 
exponential of Eq. (fTT|) . we may therefore replace the 
DDI potential as Ai2{z) -2C12I {Tiw'^) 5{z). We then 
obtain {^!i{L,L/v)) = q;i(0) exp(i(/)i2n2), with (^12 given 
by Eq. ([9]). This indicates that, at the output from the 
medium ['I'i(L + 0) = £i{L + Q)], the coherent probe field 
has acquired a phase proportional to the number of pho- 
tons 7i2 in the signal field. This phase can be measured 
by, e.g., a single-port homodyne detection using another 
coherent field of the same amplitude \ai\. The average 
detector signal is then s{n2) — 4\ai\'^ siv? [(j)i2n2 / 2) with 
the corresponding uncertainty 5s{n2) — y2s(n2)- Our 
aim is to distinguish with high probability the photon 
number states with 712 S [0,nmax]- This requires that 
0i2"Tnax < while the measurement uncertainty con- 
straint yields 5(712) — 5(712 — 1) > ^[1^5(772) + 5s{n2 — 1)]. 

Above we have neglected the self interaction within 
the probe pulse, which would otherwise dephase the 
coherent state. We can estimate its effect as fol- 
lows: - sin^ 9 / L Jl^'^dfjl^dz' Aii{z' -vt')\ai{z' -vt')\'' ~ 
2CiiSvii^ 9 / {fiw^v) |q!i(0)P, which should be small com- 
pared to 012. This leads to the condition 2Cii|q!i(0)P < 
Ci2- (Note that, as long as we are concerned with de- 
termining the photon number in the signal field, its self- 
interaction is immaterial.) Thus, for the QND measure- 
ment of the signal photon number by a coherent probe, 
the Rydberg states \d1_2) should be chosen such that 
pdi < Pd2/(2|ai(0)P), and therefore the self-interaction 
of the probe, Cn oc p\^, is small compared to the cross 
interaction C12 c>c pdipd^- On the other hand, to real- 
ize the CPHASE gate between two single photon pulses, 
012 = TT, both states |c?i.2) should have large and com- 
parable dipole moments pdi 2 so that C12 is large. 
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Of course, in all cases we need to satisfy condition (ii), 
since otherwise the DDI frequency shifts beyond the EIT 
transparency window would induce strong self and/or 
cross absorption of the fields [l^. We therefore require 
that 



2Cw sin^ e 
hw^L 



max(ir (z)) < Suji {I, I' = 1, 2). (12) 



In terms of the phase shift per photon (pw, Eq. (jH]), and 
assuming smooth n;-photon pulses of lengths Tiv < L, 
we then have 2(f>iii ni' < Ti'SoJi for cross- interaction and 
2(j)ii{ni — 1) < TiSiOi for self interaction. In turn, the 
product of the pulse duration and EIT bandwidth is re- 
stricted by the optical depth as Ti'Suji < ^\/kiL. We 
thus obtain that the maximal cross and self phase shifts 
are limited by 



4>U' nv, (t)n{ni - 1) < 



(13) 



Alternatively, the photon number in each pulse is limited 

by 



ni < 



(14) 



To relate the foregoing discussion to a realistic experi- 
ment, we assume a hollow-core waveguide of length L ^ 
1cm with the lowest transverse mode of width WfC^i 2/im 
[iMl. The waveguide is filled with iV ~ 5 x 10'' cold Rb 



atoms tightly confined by a guided dipole trap to Wa — 
2/im, leading to the effective density p ~ 2 x 10^^ cm"'^. 
For the two quantum fields tuned to the Dl and D2 tran- 
sitions \g) — 7> |ei,2) (Ai = 795 nm, A2 = 780nm), the cor- 
responding optical depths arc kiL ~ 600 and K2L ~ 580. 
With 7gei 1.8 X 10^ s-\ 7ge2 ^ 1.9 x 10^ s^^, and 
taking 0.1 ~ 7.35 x 10*^ rad/s, 0.2 ~ 7.43 x 10^ rad/s, 
the group velocities are vi^2 = 100 m/s. The bandwidth 
of the pulses Tj^^ '> v/ L — 10"* s^^ is smaller than the 
EIT bandwidth 6uJi ~ 1.2 x 10^ rad/s. To realize the 
CPHASE gate, we choose the Rydberg states |c?i.2) with 
Pdi = pd2 = 315eao (quantum numbers n = 15 and 
q — 71—1), leading to the conditional phase shift 4>i2 = tt. 
For the QND measurement of photon number ^2 < 2 in 
the signal field with a weak coherent probe ~ 4, the 
corresponding dipole moments for the Rydberg states are 
pdi = 50eao and = 450eao, leading to the cross-phase 
shift per photon of 0i2 = 0.7. We have verified that in 
both cases the DDI frequency shifts are within the EIT 
window 5uji [cf. Eq. (fT2|) ]. 

Hence, the present scheme enables a realization of de- 
terministic quantum gates with photonic qubits and is 
capable to distinguish with high probability the photon 
number states via QND measurement, which can serve 
as a heralded source of single of few photon pulses. In 
closing, we note that all the necessary ingredients of 
our proposal, including EIT via Rydberg states [H, [lB| 
and in hollow-core waveguides |17H19| , have already been 
demonstrated experimentally. 
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